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Abstract 





An analysis is presented for the flow of non-Newtonian Williamson fluid due to permeable sheet. The 
recent modification of Fourier’s law of heat conduction, i.e. Model of Cattaneo-Christov heat flux is 
studied in present article to analysis the heat flux. The permeable sheet undergoes suction and the effect 
of Soret and Dufour is also investigated The governing boundary layer equations are transformed into a 
non-dimensional form using appropriate set of transformation and the resulting non-linear coupled 
system of ordinary differential equations are solved using robust Homotopy Analysis method. The 
influence of physical parameter affecting the velocity of fluid, heat and mass transfer are discussed and 
outputs are expressed in figures and tables. Some interesting results are found from the study that the 
velocity profile is decelerating with increase in magnetic parameter. Increasing Darcy number 
accelerates the flow due to a corresponding rise in permeability of the regime and concomitant decrease 


in Darcian impedance. 
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1.Introduction 


The electrically conducting fluid flow in the 
presence of magnetic field (Magneto- 
Hydrodynamics) over permeable stretching sheet 
has importance and applications in heat exchanger, 
metal-working process in hot rolling, cooling 
process of nuclear reactor containment vessel and 
petrochemical industry. Crane [1] had provided 
analytical solution for flow of incompressible 
Newtonian fluid over a stretching sheet and the 
study of magneto-hydrodynamic (MHD) was first 
done by Alfven [2]. By incorporating MHD effect, 
the impact of boundary layer flow due to 
stretching sheet can be improved by applying 
magnetic field normal to the plane of the sheet 
which has the electrically conducting fluid. 
Further, the flow of electrical conducting 
Newtonian fluid in presence of uniform magnetic 
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field between two parallel plates is analysed by 
Hartman and Lazarus [3]. Kumaran et al. [4] 
analyzed the electrically conducting fluid flow past 
a stretching permeable sheet. Nayak et al. [5] 
investigated the natural convective due to flow of 
nanofluid over permeable sheet with thermal 
radiation in presence of magnetic field. Recently, 
Khashi’ie et al. [6] inspected the influence of 
stratification on MHD stagnation point flow past a 
stretching surface embedded in a porous medium. 

Heat transfer phenomenon in diverse practical 
applicants used to predict by the empirical formula 
of Fourier’s law [7]. However, Fourier’s model has 
certain short comes that any disturbance felt 
instantly spread by whole material. To overcome 
this short comes of Fourier’s model, Cattaneo [8] 
modified the Fourier’s model by adding thermal 
relaxation time for heat flux. Further, Christov [9] 
used Oldroyds derivative in order to preserve 
material invariance in Cattaneo’s model. Many 
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researchers have used Cattaneo-Christov model to 
analyse the heat flux. Recently, Vasu et al. [10] 
examined the natural convectional flow of Jeffrey 
nanofluid by Cattaneo-Christov model past vertical 
plate. Reddy et al. [11] have studied influence of 
Cattaneo-Christov model for hybrid dusty fluid 
flow. 

The applications of non-Newtonian flow are well 
known in science and engineering especially in 
field of bioscience (blood flow, stenosis) and 
industry like melts and solution of polymer with 
high molecular weights, emulsion coated sheets 
like photographic films, extrusion of plastic sheets 
etc. A single constitutive equation is not sufficient 
to explore the different behaviour of fluids. Many 
rheological models [12-14] have been developed 
time to time to overcome such limitation. The 
power law fluid [12], Spriggs fluid [13], Yasuda- 
Carreau fluid [12] and Powell-Eyring [14] etc. are 
some fluids which explain the properties of shear 
thinning and viscoelastic fluids. One of such fluid 
is Williamson fluid [15] which is shear thinning 
fluid. Williamson [15] proposed one of the shear 
thinning models to analyse the flow of shear 
thinning materials and _ verified the results 
experimentally. Many researchers have 
contributed valuable work to this field in past few 
years. Nadeem et al. [16] scrutinized the flow of 
Williamson fluid over stretching sheet. Kumar et 
al. [17] have analyzed the impact of non-uniform 
source of heat on the Williamson fluid flow past 
curved sheet. 

Motivated by above’ discussed numerous 
applications, the present article probes the 
boundary layer flow of Williamson fluid induced 
by permeable stretching sheet in the presence of 
electro-magnetic field, chemical reaction, Soret 
and Dufour effects. Cattaneo-Christov heat flux 
model emphases the role of relaxation time and 
studies the heat transfer phenomenon. The focus is 
on the influence of magnetic field and thermal 
relaxation time on the flow velocity and 
temperature profiles. The governing equations are 
reduced to system of non-linear coupled ordinary 
differential equation by using suitable similarity 
transformation discussed in mathematical 
formulation section. Homotopy analysis method 
(HAM) [18] is used to analyse the flow behaviour 
of Williamson fluid due to permeable stretching 
sheet in presence of magnetic field. Results 
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obtained are concluded in the Conclusion section. 


2. Mathematical Modeling 


Considering two-dimensional laminar and steady 
flow of Williamson fluid over the sheet y=0, 
provided that the flow is being restricted to y>0. 
In order to produce the stretching, equal and 
opposite forces are applied along xaxis. Origin is 
fixed. Heat transfer follows the Cattaneo-Christov 
model. The effect of Soret and Dufour is also 
considered. The equation which governs the flow 
generated by stretching sheet, that is continuity, 
momentum, energy and concentration equations, 
are given by 
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And the boundary conditions are 
At y=0, u=u,, v=v,,T=T,, C=C, (5) 
yoo, u=0, T=T,, C=C, (6) 





whereu(x, y) and v(x,y) are the velocity 
components along the flow direction and normal to 
the flow direction. Further, T is temperature, Tyis 
the temperature at the wall of the sheet and v is the 


kinematic viscosity. [is fluid parameter. v, is 
mass flux velocity. Noting that if v,>0 
corresponds to injection while v,, <Ocorresponds 


to suction. 6, and £. are coefficient of thermal and 
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concentration expansion. , a is the thermal 
diffusivity and h is the heat transfer coefficient. 
Introducing 


1/2 Uu a T -T. 
=(vxu ; =| , On)= — 
y=(vmu,) f(y). 1 [“*) yi: A= 
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wherey is the stream function defined asu = ne 

Y 
Oy ..,. ; i ae 
and v= "9 which identically satisfies continuity 
Xx 
equation.. Thus, we have 
_— df = 1/2 

u= ar and vy = —(ax)"~ f (77) (7) 


where prime denotes differentiation with respect to 
yn. After introducing (7), governing equations (1)- 
(4) and corresponding boundary conditions, (5)-(6) 
are reduced to 
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where prime denotes differentiation with respect to 
77, and material fluid parameter is represented by 
We, Pris the Prandtl number. Aare fluid 
parameters. /,is thermal relaxation parameter. D, , 





_BCy-Ca) Dkr Ty To) py _ Dkr(Cy-Cu) 
BAT, -T,)) TC, -C,)° ec, (T, -T.n) 
and S are respectively Dufour 
number,concentration buoyancy parameter, Soret 
number, Darcy and suction parameter. 
The physical quantities of interest are the skin 
friction coefficient C,and the local Nusselt 


number Nu, which are given by 
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C — Ty ; ‘ Nu, = Xq,, 


and 
Sh, =——_*—— 

D(C,, -C,,) 
After transformation, we have 
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then dimensionless shear stress and rate of heat 
transfer and rate of mass transfer becomes 
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3. Methodology 
Homotopy analysis technique [18] is semi 


numerical strategy and utilized by many scientists 
and researchers. This strategy is applied in 
numerous fluid flow problems just as applied 
science and engineering. The governing equations 
(1)-(6) are transformed to non-dimensionalized 
form (9)-(13) by using transformation (7). The 
resultant system of differential equations (9)-(13) 
is finally solved via HAM. Now to obtain the 
series solution of (9) - (13) using HAM, we choose 
the appropriate initial guesses f, ,6, and ¢ off, 
OA and ¢ in the following form 

fo =S—l+e”, 6,=e°7 and ¢ =e" 


Symbolic computation software MATHEMATICA 
9 is used for computation purpose. The HAM 
depends highly on convergence control parameter 
tt, we plot the #-curves to get the admissible 
range off, #,and #, forf, gand hrespectively 
shown in Fig. 1 and Fig. 2. Also Table 1 shows the 
comparison of present results with the previously 
published result [19] and found to be in good 
agreement. The appropriate range for the control 
parametersh,, h, and h, are—0.9<%, <—0.28, 


~0.6 < tr, <—0.25and—0.65 < #, <-0.25 at 16” 
order of approximation for values 


Le =2, « =0.55, Nr =0.25, Da =2, K, =0.2 
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Pr =7, Sr =0.25, 4, =0.1, S =0.1, D, =0.25, 
A=0.5, 7 =0.1, M =0.1, 
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Table.1.Comparison of Present Result 
~ Ty, 
Ny Present Nadeem et al. 
[19] 
0 1.005045 1 
0.1 0.97041 0.976588 
0.2 0.93734 0.939817 
0.3 0.86541 0.88272 
4. Result and Analysis 
The stream flow of an electrically leading 


Williamson fluid due to horizontal extending 
surface embedded in a porous medium has been 
considered with the impacts of the Cattaneo- 
Christov heat flux model, Soret and Dufour 
boundaries. Here an endeavour is made to separate 
an insightful answer for the exceptionally 
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confounded and profoundly coupled arrangement 
of conditions. Fig.3 depicts the variation of 
velocity with 7 for different values of Hartman 
number (magnetic parameter). It is found that as 
the effect of magnetic field increases, the velocity 
of Williamson fluid fall down significantly. It is 
due to drag force produced by magnetic field to the 
electrically leading Williamson fluid which results 
in lower down the velocity of the fluid within the 
boundary layer. 


— Williamson Fluid 


Newtonian Fluid 





0 2 n 4 6 
Fig.3.Variation of velocity for different M 


The variation of temperature profile with 7 for 
various values of MHD parameter is explored in 
Fig.4. It can be seen from the Fig. 4 that the 
thermal boundary layer increases with rise in 
MHD parameter. The influence of Lorentz force 
(drag force) reduced the velocity and generates 
friction on the flow which causes more heat 
energy. Due to this friction, temperature 
distribution in the flow field is increased. The 
behaviour of Williamson and Newtonian fluid are 
same for the effect of Hartman number on the 
profiles of velocity, temperature and concentration. 
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Fig. 4.Variation of temperature for different M 
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The effect of thermal relaxation parameter on the 
temperature field is explained in Fig. 5. A, =Ois 


corresponds to the Fourier’s model. It is observed 
from the Fig. 5 that due to relaxation parameter, 
the temperature decreases. Also temperature for 
Cattaneo-Christov heat flux model is less than that 
of Fourier’s model. 





Fig. 5 Variation of temperature for different 
thermal relaxation parameter 


Conclusions 


The presented analysis leads to the following 

major outcomes: 

e Convergence of the homotopy solution is 
discussed and the present results are compared 
with the existing results in published articles. 

e Velocity profile is a reducing with magnetic 
parameter M due to effect of Lorentz force 
whereas the temperature field is an increasing 
function of magnetic parameter M because of 
friction between the wall of the sheet and fluid 
or within fluid causes heat generation. 

e Temperature distribution is more in_ the 
Fourier’s model for heat flux as compared to 
Cattaneo-Christov heat flux model. 

e Temperature and thermal boundary layer 
thickness are decreasing functions of relaxation 
time 

e The behaviours of parameters in Cattaneo- 
Christov model are qualitatively similar to 
those in Fourier’s heat conduction law 
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